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ABSTRACT 


An Expert System based on Logic Programming for thermal 

design of Shell and Tube Vertical Condensers having single shell 
pass and single tube pass has been developed. The present system 
makes use of Expert System Shell ’VIDHI’ based on backward 
chaining. 

The design process of vertical condensers has been 
transformed into a knowledge base comprising of rules, facts, 
questions and explanations. The modified ’Bell -Delaware Method’ 
has been used for the shell-side analysis and ’Kern’s Method’ for 
the tube-side analysis. The shell-side analysis considers the 
non-ideal cross flow as it occurs in a shell and tube heat 
exchanger and the suitable correction factors for the same are 
calculated to take into account the deviation from the ideal cross 
flow. For the subcooling and the desuperheating zones, average 
heat transfer coefficients have been used, while for the 
condensation zone Stepwise Method has been used, which. takes into 
account the variation in heat transfer coefficient due to 
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1. INTRODUCTION 


1 . 1 Artificial Intelligence And Expert System : 

Artificial Intelligence (AI) investigates symbolic, 
non-algorithmic reasoning processes and the representation of 
machine inference. AI seeks to transform the way we interact with 
the computers, enlarge the ’cognitive’ abilities of the computers, 
thereby opening up many new application areas like EXPERT SYSTEMS 
etc. , besides making efforts to provide a computational theory of 
intelligence. 

An expert system is a program containing knowledge from 
a restricted domain which uses complex inferential reasoning or 
logical decision making process to carry out tasks normally 
performed by human experts. An expert system must also be able to: 

* engage in a dialogue with a user to acquire the relevant 
details of the problem 

* explain the problem solving process 

* get easily modified as new discoveries are made 

* deal with partial information. The solutions should degrade 

gracefully rather than suddenly and completely. 
(Sangal,1985) 

Typical applications of Expert Systems involve 
interpretation of large amount of data ( possibly repetitive ) or 
problems requiring many years of experience and a large store of 



unususal facts associated with that experience. Thus, an expert 

system is based upon encoding the human expert’s knowledge. 
(Engman , 1986 ) 

1 . 2 Structure of the Expert System : 

In developing an expert system, the programmer programs 
the computer on what to do without telling it how to do it. This 
is a major demarcation from historic computer programming, which 
is typically a structured set of instructions telling the computer 
how to process a set of data. Though it has historically provided 
a powerful and fast analysis tool for specific problems but not a 
very flexible tool for handling missing data, deductive processes, 
or poorly defined parameters. 

The present expert system is based on the most popular 
rule-based system or production system. It consists of three major 
components : 

* a knowledge base consisting of if-then rules 

* a current context or facts pertaining to the particular 
problem being solved by the system 

of 

* an interpreter deciding the order of application rules . 

1 . 3 Present Work : 

The objective of the present work is to build an expert 
system for the 

* selection of proper type of condenser and 

* thermal design of vertical downf low(both shell-side and 
tube-side condensation) and vertical upflow (tube-side 
condensation ) condensers . 



EXPERT USER 



PRODUCTION SYSTEM 


Fig. 1.1 : Structure of an Expert-System 
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1 . 4 Motivations Behind iiie E T . e .g ent Work '■ 

It is seen from literature available on the design of 

condensers that a very large number of decisions involving logic 
have to be made in the actual design, e.g. choice of a suitable 
type of condenser , deciding the type of the bundle, the type of 
pitch layout etc. .Many times certain parameters are already known 
and so need not be computed and overall design will depend on 

these known values . The final design is obtained by trial and 

error until certain criteria are satisfied. There . is no unique 

solution to the design problem, whichever solution suits the needs 
of the user is ultimately the chosen solution. It is for such 
ill-structured problems requiring the search in solution space 
that expert system is most suited. 

Thus the important motivations behind the present work 

are : 

(a) to increase the efficiency of an expert 

(b) to incorporate additional information in this area . 

(c) to provide an aid to expert for confirming his decisions 

(d) to act as teaching aid. 

1.- 5 , Organisation si. the The g l ^ : 

Chapter 2 discusses the Expert System Shell ’VIDHI’ 
which is the basis for the present expert system. Chapters 3 and 4 
give all the relevant details for the selection of the type of 
condensers and the thermal design of vertical condensers. Chapter 
5 discusses the structuring of the knowledge base used in the 
expert system developed. Chapter 6 illustrates the use of the 
system and some sample sessions have been recorded. Limitations of 
the system and scope for further work is also outlined. 
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2. EXPERT SYSTEM SHELL ’VIDHI 


2 . 1 Introduction : 

An expert system shell provides the user with a 
provision for efficient storage and retrieval of assertions, an 
appropriate user interface, ability to pose questions to the user 
(when facts and rules fail) etc., (Sangal, 1988) . Knowledge related 
to any particular domain can be entered into it to develop the 
expert system in that domain. 

The shell 'VIDHI’ is written in LISP .So while using it 
one can make use of all the features and facilities of LISP, 
However, use of VIDHI provides a simpler syntax than that which 
would be necssary if the program were to be written in LISP. 

2 . 2 Pattern-Matching : 

A pattern is an S-expr which might contain zero or more 
occurrences of wild cards starting with the character ’?’. A data 
item is an S-expr which has no occurrences of wild cards. A 
pattern matches a data item if on replacing each of the wild cards 
by an appropriate S-expr, we get the data item.For example, the 
following two matches: 

Pattern : (fluid-name hot-fluid ?name-h) 

Data-item: (fluid-name hot -fluid steam) 

because on substituting ?name-h in the pattern by 
’steam’ we get the data item. 



2 . 3 Logic Programming : 


Logic programming discribed here is used for solving 
problems involving objects and relationships. To express 
relationships between objects, predicates are used. For example, 
to express that the pressure of condensing ammonia is high, we 
make use of two place predicate ’pressure’: 

(pressure ammonia high) 

Here ’pressure’ is the predicate and ’ammonia’ and ’high’ are 
its arguments. 

2.3.1 Terms and Formulas: 

Terms occur as arguments of the predicates in formulas 
and are either variables or constants or function-argument 
combinations, e.g. 

(1) ?x, ?person, ?p3 

(2) 3, 5.6, HIGH 

(3) a list of the form 

(<f> <pl> <p2> <pn>) 

where <f> is a function symbol (a constant) and <pl> to 
<pn> are argtiments (terms), e.g. 

(add 2.0 (product 2.0 4.0)) 

where 'add’ is a function followed by its two arguments. 

A formula can take one of the following forms; 

(1) An atomic formula is a predicate-arguments combination where 
the predicate is a symbolic atom and the arguments are terms. For 
example, 

(density water 1000.0) 

is an atomic formula having ’density’ as a predicate and, ’water’ 
and ’1000.0’ as its arguments. 
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(2) A horn-clause (formula) is of the form 

A n >- 0 <2.1) 

B <- A1 An n >- 

where B and A1 to An are atomic formulas. B is called 
the and A1 to An the If the antecedent is 

empty, (n=0) , horn clause reduces to an atomic formula. 

2.3.2 Assertions, Rules, Facts and Query- 

Formulas described in the previous section take truth 

values. Formulas known to be true are known as assertions, e.g. 

(greater 3 0) 

(= 3 (add 2 1)) 

If a horn-clause of the form given in Eq. (2-1) 
included in the data-base, then B is true whenever each member of 
A1 to An is true. The consequent formula is unconditionally true 

if “the an'teceden't is emp"ty > e.g. 

(cond-type tube-side-condensation) <" (condensing-f luid 

corrosive) 

Here (cond-type tube-side-condensation) is true xf 

(condensing-f luid corrosive) is true. 

An assertion with antecedent ( s ) is called a rule, while 

an assertion with an empty antecedent is called a fact. 

A formula whose truth value is to be determined is called 

a query or a goal formula. It is written as 

goal (cond-type ?1) 

If the goal formula has a variable, then the binding of 
a variable for which the formula true, is the answer. 

Inference- 

There are two rules of inference using which 


2.3.3 



value of a formula can be determined to give an initial set of 
assertions : 

a. Universal Instantiation ; 

Substituting a variable in a rule yields an assertion. 

b. Modus Poens (Application of a rule) : 

If the formulas in the antecedent of the assertion are 
true, then the consequent is also true. For example, from the 
following: 

(fluid-name ?x) <- (condensing-f luid ?x) 

(condensing-f luid steam) 
it can be inferred that 

(fluid-name steam) 

is true. 

2.4 Present Shell ’ VIDHI' at the User Level : 

VIDHI is based on backward chaining where one begins 
with the goal and goes backward towards the facts. Computer 
programming in VIDHI consists of 

declaring some facts about objects and their relationships 

* defining some rules about objects and their relationships 

* asking questions about objects and their relationships 

The way of using the features and characteristics of 
the ’Shell’ are described below in brief: 

2.4.1 Facts and Rules: 

The syntax for a fact is 

(def asrt <f act-name> ( <predicate> <argl><arg2> . . . <argn> ) ) 
and that for a rule is 

(def asrt <rule-name> ( <predicate><argl><arg2> . . . <argn> ) <- 


8 



( <predicate><argll ><argl2> . . .<argln>) 
( <predicate><arg21><arg22> . . . <arg2n> ) 


( <predicate><argml><argm2> . ..<arginn>)) 

2. 4. 2. Posing Questions: 

This is a mechanism in the ’Shell’ used for eliciting 
the information from the user. This is accomplished by introducing 
rules that cause questions to be issued. Such rules make use of a 
built-in predicate called ’ ask-user’ . The database is not searched 
to test the truth value of a formula involving ask-user, instead 
the interaction takes place with the user. The syntax for 
’ask-user’ is 

(ask-user (source <varl>. <var2> . . . <varn> ) 

(target <varl> <var2> ...<varm>) 

(question <text of the question>) 

(types <specl> <spec2> . . . <specm> ) ) 

All the source variables should have variable free 
values, but none of the target variable should have a value. The 
further details about ask-user can be obtained from VIDHI-Manual . 

The results obtained using ’ask-user’ are stored 
automatically. 

2.4.3 Storing Intermediate Results: 

After a goal or a subgoal has been inferred, it is 
stored in the data-base if the predicate used in that goal or 
subgoal has been declared ’intermediate’ using the following 
syntax:: 

( def intermediate <predicatel ><predicate2> . . . <predicate n> ) . 
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2.4.4 Elaborat-ion '• 

Text may be associated with a predicate providing a 
helpful explanation to the user of expert system. It is retrieved 
gand displayed when user types ’what’ in response to a query. The 
syntax i s : 

( def elaboation <predicate> ( <elaboration/explanation> )) 

2.4.5 Computational Predicates and Functions : 

Computational predicates are ’computed’ (i.e., applied 
to their arguments) to determine the truth value of a formula in 


which they occur. 

All 

the 

LISP 

computational 

functions and 

predicates can be 

used 

in 

VIDHI 

after making 

the following 


declarations: 

(defcompfunc <VIDHI-name> <LISP name>) 

(def comp <VIDHI-name) <LISP-name>) 

The use of computational predicates and functions places 
a restriction on the choice of input and output arguments. All the 
argtiments of a computational function must be known. 
Moreover, the antecedents become procedural. 

2.4.6 Functional Dependencies: 

Many times the value of an argtnnent of a predicate can 
be functionally determined by the value of other arguments. For 
example, density of water is unique at a given temperature . Thus 
in 

(density ?temp ?dens) 
the ?temp determines ?dens and 

(density 100.0 958.3) 

(density 100.0 1000.0) 



should not be allowed to be present simultaneously. 


This idea is captured using declaration of the type : 
(deffd <prdeicate 1> <prdeicate 2> . . . <predicate n>) 

It is assumed in the current implementation that only the first 
argument is independent and all other arguments are dependent on 
the first one. The declaration of functional dependency allows 
inconsistencies to be detected. 

2.4.7 ’retract-all’ Predicate: 

Many parameters in any design problem have some fixed 
value and once this value is determined, the predicates associated 
with these parameters are no longer needed. In fact there presence 
may cause problems during backtracking. So the rules involving 
these predcates are made inapplicable whenever the following 
built-in-predicate appears in the antecedent of a rule: 

(retract-ali <predicate>) 

2.4.8 Inference Mechanism: 

The inference mechanism of the system is similar to that 

I 

of PROLOG. It doesn’t distinguish between rules, facts or 
question-rules (the rules involving an antecedent having ask-user 
predicate). 

The system accepts any query when the function ’goal’ is 
used, e.g. if the query is (area-req ?1) then 

’goal(area-req ?1)’ 
will try to satisfy the query. 

The infer algorithm of the ’Shell' takes the query and 
does a series of operations : 
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1 . It searches the data-base and checks whether the predicate 
currently involved is a computational function or predicate and if 
so, then the arguments are passed onto the compute predicate or 
function after substituting the values for the arguments of the 
predicate . 

2. If the predicate is one of ’not’, ’ask-user’, or ’retract-all’ 
then the actions are taken as described in the previous section. 

3. If the predicate is not of the type (1) or (2) then the 
data-base is searched for the rules and facts and stored for it. 
these are tried in the order in which they appear in the 
data-base. For facts no new subqueries are generated while for 
rules a new subquery is generated for every predicate in the 
antecedent . 

The value of arguments of any predicate are stored 
if the predicate (that gets satisfied used} in that goal or subgoal 
has been declared ’intermediate’. 
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3. PRINCIPLES AND CLASSIFICATION OF CONDENSERS 


3 . 1 Modes of Condensation : 

Condensation occurs when the temperature of vapour is 
reduced below its saturation temperature. In industrial 
equipments, the process commonly results from contact between the 
vapour and a cool surface and is called ’Surface 
Condensation’ . There are two forms of Surface Condensation (Fig. 
3.1) : 

(i) Filmwise condensation 

(ii) Dropwise condensation . 

(i) Filmwise Condensation : 

The condensate forms a continuous film on the cooled 
surface. This is the most important mode of condensation occurring 
in industrial equipments and is discusssed in detail in this 
Chapter . 

(ii) Dropwise Condensation : 

IS 

This occurs when the condensate Aformed as droplets on a 
cooled surface instead of a continuous film . High heat transfer 
coefficients can be obtained with dropwise condensation, but this 
is difficult to achieve and maintain in heat exchangers. 

3.2 Resistances to Condensation : 

Figure 3.2 shows the filmwise condensation of pure 
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(a) (b) 


Fig. 3.1 Modes of Condensation 

(a) Film Condensation 

(b) Dropwise Condensation 



Wall Condensate 
film 


Fig. 3.2 Resistances to Condensation for the Case of a Pure 
Vapour 
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vapour which has been assumed as saturated. Temperat-ure drops exist 
across the condensate film , through the wall and from the wall 
into the coolant (Butterv.orth, 1935) . 

Heat transfer coefficient h for the condensate film can 

C 

be written as 


h = 
c 




T -T 
^Sat 


(3.1) 


where 


heat flux through the film 


Tj = interface temperature 


'Sat 


saturation temperature of vapour 


= wall temperature 

In general the following parameters affect the film 
resistances - 

(a) Geometry of surface: (i) Inside of horizontal or vertical tubes 

(ii) Outside of horizontal or vertical 
baffled tube bundles 

(b) Vapour hydrodynamics :( i ) Condensate film laminar (no vapour 

shear) 

(ii) Gravity induced turbulency 

(iii) Vapour shear induced turbulency 
(i)Single component 

(ii)Multi- component fully condensable 
mixtures 

(iii) Vapour with noncondensable gas present 
The Nusselt number for condensation may be expressed as 


(c) Substance: 


Nu = 
o 


h 6 


k-, 


X 10 


where , h = heat transfer coefficient for condensation 
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= thermal conductivity (liquid) 

6 = film thickness ; it depends on the 

force of gravity & vapour shear force 


3 . 3 Condensation Inside or Outside Vertical Tubes ; 

Value of Nusselt number in this case will depend upon 
the value of Reynolds number and Prandtl number where , 


4r. 


Re^ = 


(3.2) 




Pr= 

l 


Si^i 


(3.3) 




M 


where 


X condensation drainage length per 


unit unit v/idth and 

S = riD, D being the internal tube diameter 
for inside tube condensation and outside 
tube diameter for condensation outside 
the tubes . 

Heat transfer values for both the cases , viz. , shear 
dominated and gravity dominated flows , should be calculated 
and maximum of these should be taken as the actual value. 

I 

3.3.1 Gravity Dominated Flow: 

For Re < 30 , the following relation can be used : 


hl= h 


Re 


-1/3 


Lp^(Pl-Pg) g^ 


(3.4) 


For 30< Re <1600 , laminar waves start building up and 
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then , 


h^=0.756 




- 0.22 


(Re,) 


(3.5) 


For Re > 1600 , turbulence starts having significant 
effect on heat transfer and then the value of h^ depends on the 
Prandtl number also . 

For Pr < 10.0 , 


h^=0.023 Ref-^^Pr®-^ 


Pi (P^- Pg) 


(3.6) 


For Pr > 10.0 , the above equation can be used with Pr = 10.0. 


3.3.2 Shear Dominated Flow: 

Henstock and Hanratty (1976)has suggested the following 
formulation for evaluation of hg- 


mg (1-y) 


Pg(nDV4) 


X 10' 


(3.7) 


where D is the tube inside diameter for flow inside tubes or the 
hydraulic mean diameter for downward vapour flow outside the tubes 
in a bundle , 


r=r(0.707 Ref 0379 Ref'®)^*^j 




f = 0.046 Re 

g g 
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G= 


p g D 
L n 


■ 2 

P V f 
g g g 


X 10 


-3 


r r) 

F = X 

D o. £> rj 

Re g 








1 + 1400 F -|l - exp 




_ (1+1400F) 
13.2 F G 


3/2^ 


} 


Ti= 0.5 f, V^‘ 


Tj = 


Pi Tj 




2/3 


Value of fhe critical Reynolds number (Re ) depends upon 


the value of Tj. 



For 

Tj > 9.04 

. Re^ 

= 50 

& 

For 

Tj < 9.04 

, Re^= 

1600 

Re 

< Re 

then , 




.+^3 


h* = 1,41 

and if Re > Re then , 
c 


Hh Hh” 

h"^ = (Tp 



r 1.41 ]\ 1 

r 0.071 Pr?*^l® 



1 ReK24 J J 


1/m 


where m - 0.5{Pr+3) 

hg is related to h^ by the following expression 


hg = h"^ k 


1 ^ . 
L Pi' 




,2 


^ ^1-^g) « 


-] 


-1/3 


(3.8) 



h can be 


Afber knovjing the values of & hg 


determined by 

he= MAX [hj^,h2] 


(3.9) 


3.4 Flooding Phenomenon : 

This phenomenon becomes important for upward vapour 
flow. For low vapour flow rates the condensate runs freely from 
the bottom of the tube . However ,if the vapour velocity is 
gradually increased, at a certain point large waves and 
disturbances occur at the bottom of the tube with the liquid being 
held up periodically by the inflow of vapour . Some liquid is 
discharged from the top of the tube , this phenomenon is known as 
flooding and the lowest vapour velocity at which it occurs is 
known as the flooding velocity . 

Flooding mass velocity is defined as , 

my= density x condensate flow velocity . (3.10) 

Holmes (1950) data can be expressed by the following 


equation - 



80 (j^000 O') 

^0.14 (COS ©)^-^^(Wj^/Wg 


.09 ^0.5 

Prr 


J0.07 


(3.11) 


Here D - tube ID (mm) 

3 

p, -densities (kg/m ) 

<y -surface tension (N/m) 

rj -viscosity (cP) 

B -taper of tube end from horizontal 

w, /w - weight ratio of condensate to vapour at 

1 g 

the bott om of the tube^^^ ^ ^ ^ 
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This flooding velocity should not exceed the above 
value of . 

3 . 5 Classification of C ondensers : 

Condensers can be classified in various different ways: 
Type of ■ Industry : Chemical process industry , Power industry 
Type of Condensation : Dropwise condensation , Filmwise 

condensation, direct contact condensation 
Extent of condensation : Desuperheater condenser. Condenser, 

Condenser-subcooler , 
Desuperheater-condenser-subcooler 
Type of construction ; Shell-and-tube, plate, spiral, 

air-cooled . 

Orientation of condenser : Vertical, horizontal 
Occurence of condensation : Inside tubes, outside tubes 
Condensation controlling mechanism : Gravity controlled 

Shear controlled 

The details regarding various types of shell and tubes 
condensers can be obtained from AgrawaK 1988) . 

3 . 6 Selection of Condenser : 

Selection of condenser types involves consideration of 
numerous conflicting requirements. The various classes of 
condensation are indicated in Flow-chart shown in Fig. 3.3 where a 
preliminary selection of a condenser type is to be evaluated. The 
Tables 3. 1 to 3.4 then give the opinion on the applicability and 
predictability of the various types of condensers. The examination 
of Flow-chart and Tables should lead to proper selection of 
condenser type . 
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£'xs> 3,3 Flow Ch.siT't for "tlio Ssloction of s Concisiissr 
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Table 3.1: Tobal Tube-side Condensation 



Horizantal 

Vertical -down 

Vertical-up 




flow 

flow 

Single component 
Vapour 

e 

coi 

^ a%j 

/ b 

Multi component 
vapour 

f 

b 

€ t> 

/ 

Sub-cooled 





Condenaate 

p 

X 

^ CLX) 

X X 

Pressure drop 





High 

e 

b 

s c 

X c 

Lov 

p 

b 

/ c 


Coolant 





t-iquid 

e 


e 

s 

OOLB 

e 


e 

e 

Boiling 

e 


€ 

e 


Comments; Horizontal condensers have possibility of slugging. 

Vertical down flow handles dirty or polymerising 
vapours. . 









Horizontal 


Vertical 
down flow 


Vertical up 
flow 


Singly component 
vapour 



Sub cooltfrd 
condensata 



Singla component 
vapour 


Multi componanl 
Vapour 


Horizantal Horisontal Vertical Vertical 
cross baffled Down Flow up flow 




e <= 


e 


p c 


Sub coolad 
Condaniaata 


/ fc 


X X 


Pressure Drop 


X X 


€ 



Gomments: Horizontal Cross can have possible venting problems 
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Table 3.4: faxiisl Shell-side CopdefiggtiQB 


Horizontal Horizontal Vertical Vertical 

cross baffled Down Flow up flow 


SiLn^Io compon^-ni 

P 

b 

p 

b 

p h 

X X 

vapour 







Multi componc-nl 
Vapour 

/ 

b 

P 

c 

X 5 

/ 

sub coolod 
Condortsaio 

e 

h 

6 

b 

X a%> 

/ 

Pressure Drop 







High 

e 

b 

€ 

b 

X c 

X X 

L,ov 

e 

b 

P 

b 

/ 


Coolant 







Liquid 

€ 



6 

e 

€ 

aoLSi 

6 


e 

e 

X 

Boiling 

X 



X 

e 

X 

Comment: Horizontal 

cross 

can 

have possible venting problems. 


^^acceptabili g - good , / Fair , p - poor 

X = not acceptahle or not recommended.. 

predictability: av = average ^ f 5 2 

= fair < % 

^ = poor ^ 50 % 

X = no method or not recommended 




4. METHCaX>LOGY FOR THE DESIGN OF CONDENSERS 


4 . 1 Introduction ; 

Overall design methodology remains the same as that for 
a STHE involving no phase change. The details regarding the design 
method of STHE without phase change can be found in Bhaskare 
(1986). Tube-side flow is well-defined and so can be easily 
computed. For shell-side Bell-Delaware-Method has been suitably 
modified to take into account the phase change part and for the 
final rating of the condensation zone successive summation method 
or Stepwise -Method has been used. The additional complications 
which arise in a condenser are as follows : 

1 . Sensible heat-transfer coefficients are rather low compared to 
condensation coefficients. 

2. The non-linear heat release requires piecewise calculations 
instead of an overall calculations as is generally done for no 
phase change units. 

3. The change in heat-transfer coefficient takes place at each 
step due to changes in physical properties of the mixture 
undergoing compositional changes (Gupta, 1986) . 

4.2 - 

The basic logical structure of the design process has 
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Fig. 4.1 Logic of the Design Method 









been shown in Fig. 4.1. It involves rating both in terms of heat 
transfer as well as pressure drop calculations ,but in the present 
thesis the rules for the evaluation of pressure drop for the two 
phase flow are not included (Bell, 1986a) 

4 . 3 Approximate Sizing of the Condensers : 

An assumed or approximate design of condenser is 
necessary to start the iteration because the values of various 
heat transfer coefficients involved can not be determined unless 
the geometrical features of the condenser are known. The steps 
involved in obtaining an approximate design are as follows : 


4.3.1 Selection of Condenser (See Chapter 3) 


4.3.2 Determination of Heat-rate : 

The heat duty which a heat exchanger is 
required to perform is given by : 


(4.1) 


4.3.3 Overall LMTD : 

By knowing or selecting the proper coolant temperatures 
the overall log mean temperature difference (LMTD) can be 
determined as follows : 

For CO -current flow : 

/ m ■ fn ■ ■ 

(4.2a) 


AT - “ ^\o" ^co^ 

, ■ 


In 


f Kj ci j 

vp T 
ho 
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For counter flow. 


AT =- 


( T, -T )- (T^ - T , ) 

Hi CO h O' C i 


In- 


T, T 
ni CO 

- T . 

Ko Cl 


(4.2b) 


4.3.4 Overall Heat Transfer Coefficient: 

This should be estimated from previous experience of the 
designer or from Table 4.1 . 


4.3.5 Approximate Area : 

The approximate area is calculated using the following 

equation : 


q 


4.3.6 Geometrical Features : 

After determining the area needed to satisfy heat 
transfer duty requirements, the geometrical features of the heat 
exchanger like tubesize, tube-length, shell-size, pitch and 
arrangement of baffles can be determined as described in the 
following section. 

4.4 Shell-side and Tube- side Param eteys. : 

The constructional features of STHE are shown in Fig. 
4.2. The surface area can be written as ; 

A = A (D , L , D. , type of pitch-layout) 


(4.4) 




!. Shell 
Z Tubes 
J. Chonnel 


10. Floating -head clomp 

11. Flooting tubesheet skirt 

12. Flonges 


4. Flooting-heod cover 13. Stuffing-bo* flonge 

5. Flooting tubesheet W- follower 

a Stotionory tubesheet 15. Supports 


I Boffles 
a Spocers 
a Tie rods 


16. Tubeside nozzles 

17. Shellside nozzles 

18. Impingement plote 


Fig. 4,2 Homenoiature for Shell and Tube Heat Exchanger 
Componen'ts 
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lab l.e 4.1 ; values ot Overall Coefficients 


vapour 

coolant 

li W / s 'Cl — m. 

alcohol 

WCLt^T 

50-1100 

dow therm 

tal 1 ail 

40-450 

dow therm 

daw t harm 

450-680 

high boiling" 
hydrocarbons - 
under vac c um^ 

watar 

100-280 

1 ow bo Hi ng 
hydrocarbons 

xjxitar 

450-1140 

organic solvents 

xiKztar 

550-1140 

hydroc arbons 

oi 1 

140-230 

herosene 

watar 

170-370 

Kerosene 

oi 1 

110-170 

naphtha 

wxter 

280-430 

naphtha 

oi 1 

110-170 

steam 

fo&d watar 

2200-5700 

vege t ab le oils 

wat&T 

110-280 

organic steam, 
azeotrope 

wator 

220-450 

air coolers 



steam 

air 

730-800 

ammonia 

air 

550-680 

light hydrocarbon 

air 

450-540 

light naphtha 

air 

400-450 

freon 

air 

340-450 


heavy naphtha 


air 


340-400 


















The type of pitch-layout is characterized by an angle 
^ shown in Fig. 4.3. The relation between the heat-transfer area 
A^, and other parameters of Eq. 4.4 is given by 

A„ = 10'“ A* [ (4,6) 

The tube-layout density parameter, A is defined as 

D 

A* = 0.78 n f-~l mm~" (4.6) 

^ (L f 

Ci= 1.0 for = 45° or 90'' 

1 tp 

C,= 0.866 for - 30°. 

1 tp 


4.4.1 Shell-side Parameters : 

The detailed description regarding bundle-types, 
bundle-shell-clearances and baffle geometry are available in the 
Refs. Bhaskare(1986) and HEDH( 1983, Section 3.3.5 and 3.3.10 ). 

To compute the coefficients for heat transfer the 
calculation of various flow areas and correction factors is 
essential. The cross-flow area at shell-centerline within one 
baffle-spacing is given by 



D 


ctl 


^p.*ff 


] mm^ 


(4.T) 


where L. = effective tube pitch 

tp, eff 

=L for 30° and 90" layout 

tp 

=0.707 X L for 45" layout 

, ■ tp ■ 

The details regarding various leakage- streams and 
Bell-Delaware Method as suggested by Tinker have been discussed in 
Bell(1986^. The various leakage areas, viz., baffle window flow 



area, bundle to shell bypass area, shell to baffle leakage area, 
tube to baffle hole leakage area and associated correction factors 
j (baffle-window correction factor), j (baffle leakage correction 

'O ^ 

factor), j^Cbaffle bypass correction factor), j^(adverse tem- 
perature gradient correction factor) and j (heat transfer 

S 

correction factor for unequal end baffle spacing)can be calculated 
by the procedures given in Bhaskare( 1986 ) and Taborek (1983). 

The shell-side dimensionless parameters are : 

1. Shell-side Reynolds Number : 

The Shell-side mass velocity is given by 

M _ 

m^ = — X 10“ kg/m sec 

m 

The shell-side Reynolds Number is given by 
Re =D m / 

m t 17 

m 

2. Shell-side Prandtl Number is given by 

Pj. 10‘® (4.10) 

* k 

4.4.2 Tube-side Parameters : 

1. Tube-side Passes : 

The number of tube-side passes is determined on the 
basis of flow velocity in the tubes. The present expert system can 
design a heat exchanger with only one tube side pass. 

Mass-velocity in tubes is given by 


(4.8) 


(4.9) 
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(4.11) 



M^x 10*’ 


Tube”side Reynolds number is given by 
Re. 


CD^-2t^) m^ 




Tube-side Prandtl number is given by 


(4.12) 


c n 

pit -3 

Pr^= X 10 


(4.13) 


4.5 Shell-side H. e . a , t Transfer Coefficient : 

The basic equation for calculating the effective 
shell-side heat-transfer coefficient is given by 


i d#al 




(4.14) 


where h. . , is the heat transfer coefficient for an ideal tube 


bank and is determined as described below. 


4.5.1 Condensation Coefficient : 

If shell-side condensation is going on, then is 

determined by the expressions given in Section 3.3. 

4.5.2 No Phase-change Coefficient : 

If tube- side condensation is going on and coolant is 
flowing on the shell-side, then in the condensation zone coolant 
(shell-side) heat transfer coefficient is obtained from the 
j-factor curves given in H1DH( 1983). The ideal tube bank factor is 


written as 



3. = J. 


K' -ir] 


(4.15) 


1.25, the fitted equations for j. are : 


For = 30 staggered : 

log^^ = 0.053 - 0.7965 x + 0.1999 

For = 45° staggered : 

log,„ 3. = 0.0676 x' - 0.8282 x + 0.3 

jlC/ I# 

For 0 = 90° inline ; 

tp 

log^^ = 0.05386 x^ - 0.7543 x 
where x = log^^ (Re^) 
h .^ , = j C m (Pr (<P )"■ 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


where {<^'^) = viscosity correction factor that accounts for the 

viscosity gradient at the tube-wall temperature versus the 

viscosity at the average bulk fluid temperature rj . The term {<P 

8 8 

is computed as follows ; 


For liquids : 


^ 0.14 


RJ 


(4. 20) 


Tube-wall temperature is taken equal to condensing 
temperature in the condensing zone. 

Viscosity of fluid at different temperatures is 
estimated using 


Tj = a T 


(4.21) 
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For gases : 


Gas 

being cooled ; 

{4> f - 1.0 

(4.22) 



% r 

+ 273.15.0-25 

Gas 

being heated : 

- 1 T + 

V 

273.15 J (4.23) 

In 

desuperheating 

zone whether 

it is shell -side 

condensation 

or tube side 

condensation. 

the shell-side heat 


transfer coefficient is calculated using Eqs. (4.16) to (4.19). 
However, If the condensation is taking place on the shell side the 
shell-side heat transfer coefficient in the subcooling zone is 
calculated using the following expressions : 

In the laminar region : 

2 2 

, g P. .. f rr’ ..1/3 1 1./3 

K=0,78 1^{ (4,24) 


where L is heat transfer length in m in subcooling zone. 
The transition region is where Rej^ > Re^ and Re^ is 
critical value of Re^^ 

Re = 2460 (Pr)"®'*^ 

U 


Here , 

2 

f H ^ 1/3 

® _ O- 2 ^ O. 34 

— j Rej^ Pr^ .... 

In the turbulent region, ReL> 2300.0, 
h. = 5.7 X 10~® ^ j 


- (4.25) 

(4.26) 


4.6 Tube -side Heat f idieat ; 


4.6.1 Condensation Coefficient : 

If it is a case of tube side condensation , then 




the 


(4.27) 



and (condensation coefficient.) is obtained using the 
expressions given in Section 3.3 . 

4.6.2 No Phase Change Coefficient : 

The heat transfer coefficient for all cases of no phase 
change except subcooling is determined by j-factor curves as given 
in Section 4.5.2 .However, the fitted equation in this case is: 

log^^j^ = -0.2624 + 3.273 x - 7.412 (4.28) 

where x = log Re. 

lO i 

The actual heat transfer coefficient is then determined 
by 

h = j (Pr [<p f (4.29) 

where C<^)*’ is to be determined for the tube-side 
fluid in a manner similar to Section 4.5.2 . 

The tube side heat transfer coefficient for the tube 
side condensation case in subcooling zone is again computed by the 
Eqs. (4.24) to (4.26) . 

4.7 Overall H&al T xflhgf^ r Coeff Id ight = 

After having computed the Shell-side and the 
Tube-side heat transfer coefficient the overall heat transfer 
coefficient can be determined by using the equation : 


r 

) 3 

r. 

'■% . ■ 

(4.30) 


This will be an average value if 
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0 = 
o 


h 


+ R 


t^o- 


r. ) 


r * 0 

o 


“f s 


r + 

o 




) + [ 




h. 


all the other 



coefficient:- t.aken as average value, e.g., in desuperheating 

and subcoC'iinK .tones. For the condensation zone, it will be 
calculated at e-,oh step since the condensation coefficient is a 
strong function of vapour fraction. 


4.0 Msrthpj fax Zmm • 

Figure 4.4 shows a condenser where condensation zone 
starts at section kA .The section M is a known section where the 
temperature of coolant, temperature of condensing fluid (= 
condensing temperature) and physical properties are known. 

At the beginning of condensation, i.e., at Section AA : 

Vapour Fraction, y = 1.0 

Consider the full condensation zone to be divided into 
steps. Then for each step : 

aon« ^outlot zon# 

6 y = ;; (4.31) 

- 1 ” 0-0 > _1 
" n “ n 


For ith step, the coolant inlet temperature T^.^^ and 

shown in the Fig. 4.4. 

(1) Average condensate flow rate per unit width for ith step is — 

C1.0 - (y + <5y/2)3 X m 


outlet temperature T . are 

c#out 


nB N 

t t 


(2) Reynolds Number for ith step is 


Re. 


4 r. 

K 


(4.32) 


(4.33) 


Knowing Re. condensing heat transfer coefficient and 
overall heat transfer coefficient can be determined. 

(3) Heat load per step is given by 


SO 



i = n+1 A 


i = i B 



A {y»1.0) i=| B- 


( = n + 1 

B {y=0.0) A- 


f 


-|-B {y»0.0) 


t 


■A (y = I.O) 


Coolont Condensing fluid Coolant Condensing fluid 


^c.out * 


Tc.{n-Tc,i 



T:,out=Tc,i 


i + l 


yj^y ,Tc,in*Tc,|+i 



yi=y 


i +1 


Coolant Condensing fluid 
Counter -current condenser 


Coolant Condensing fluid 
Co- current condenser 


(b) 


Fig. 4.4.' Stepwise method for the condensation zone 
0 . Schematic of the condensation zone 

b. Representation of i^^ step indicating coolant temperatures 
and condensing -fluid’s vapour fractions 




(4) At the inlet of ith step is known. At the outlet of 


step •• 

Aq 

T “ T - 

c,oyt c,in Ijj c 

c pc 

or, for ith step— 


(4.35) 
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pc 

where negative sign is for co-current flow and positive 
sign is for counter flow. 

(5) Knowing Aq^ from step 3, area for any element can be 
calculated : 


A q„ 

AA, = -y — 

1. CM, L 


(4.37) 


where AT, is determined from Eqs. (4.2a) and 


(4.2b) . 

(6) After calculating AA^ for each step, the area required for 
condensation zone is given by : 

n 

^ (4.38) 

con mm ^ X ^ 

i S! 1 


4.9 Overview of Rating Erafisdars C al culat i ng Ars& Required : 

The rating problem can be expressed as in Table 4.2. 
Rating, in general, involves following steps : 


4.9.1 


Besuperheating Area Calculation : 

Desuperheating heat load is given by : 
q. = m c T ) 

Desuperheating zone will be at coolant inlet 


(4.39) 
point for 



Table 4.2 ■- Rating Problem 


Rating Problem 

Given : Process Specifications 
Flow rates 
Inlet temperatures 
Outlet temperatures 
Physical properties 
Fouling characteristics 
Geometry of Exchanger 

Shell inside diameter 
Outer tube limit 
Tube diameter and layout 
Baffle spacing and cut 
Calculate : Area required 



co-current heat 


exchanger and coolant outlet point in 

counter-current heat exchanger. Fig. 4.5 shows the schematics of 
co-current and counter-current heat exchangers. The terms 
co-current and counter-current just indicate the relative 
positions of the inlet and the outlet points of both the fluid 
streams and not the parallel flow. T . ,T , ,T and T 

ct con,i con^o co 

have been indicated in the Fig. 4.5 . 



LMTB correction factor for one tube pass and one shell 
pass heat exchanger is 1-0 .The detailed description of it is 
given in Bhaskare ( 1986 ) . How, desuperheating area is determined 
by : 




Condensing fluid Coolant Condensing fluid Coolant 


(a) 



Length of Condenser — ► Length of Condenser — — 

(b) 

Co- current Condenser Counter- current Condenser 


Fig. 4.5(a) Schematics of Various Heat Transfer Zones 

(b) Temperature Profiles for Co -current and Counter - 
current Condensers in the absence of Non- 
condensables. 
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(4.42) 




"dosup " u . AT F 

a#fiup JLMd& ^ 


where will be determined based on the procedure 

given in Section 4.7 . 


4.9.;? Condensation Area : 

The condensation area is calculated by the 
stepwise method as given in Section 4.6 
4.9.3 Subcooling Area : 

The subcooling heat load is given by : 

%uh ~ (4.43) 


The temperatures and so the LMTD can be 
determined on similar lines as for desuperheating area in Section 
4.7.1 . 


The overall heat transfer coefficient U ^ , for 

eubcool 

subcooling zone can be detrmined based on the procedure given in 
Section 4.5 .So the subcooling area is given by : 


^muh 


«'ub 


0' , AT , F 

«ybcoo4 JLMsub c 


(4.44) 


4.9.4 Total Area Eequired : 

The total area required from the rating program is given 
hy *. : • 

A = A, + A + A . (4.45) 

d#iiup , ■con sub 



5- RUL.K SETS OF THE EXPERT SYSTEM 


5 . 1 FJ-S'W f-QX Ihe • 

Fig. 5.1 gives the flow-chart in brief for the design of 
condensers. The steps to be followed can be summarised as given 
bolow : 

1 . Decide about the type of condenser. This involves choice 
between shell -side and tube-side condensation and the selection 
regarding vertical or horizontal, co-current or counter- current 
(both pertaining to E-type shell) heat exchanger. 

2- Estimate the approximate area from the property-data and 
specifications . 

3. Decide about the type of bundle, baffle geometry and shell and 
tube geometry from the approximate area and the limitations 
imposed upon the design. 

4. If the choice is vertical-upf low condenser involving tube-side 

condensation, then ensure that the tube-side mass velocity is less 

than the flooding mass velocity else modify design by going 

back to Step 3 . 

5 . Determine the various correction factors from the geometry of 
STHE and various flow areas. 
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6 . Now calculate desuperheating area after estimating the overall 
heat transfer coefficient for desuperheating zone. Area required 
for the condensation zone should be calculated after estimating 
the overall heat transfer coefficient for condensation for each 
step (stepwise method). Area required for subcooling zone is 
calculated in the similar manner as for the de.superheating zone. 

7. Total area required is known from Step 6 . Compare this actual 
required area with the area available for heat transfer. If the 
comparison is acceptable then one proceeds with further design; 
else the design should be modified by going back to Step 3. 


The 

rule 

set 

for 

the expert 

system 

is, however, 

developed by 

using 

the 

tree 

structure 

for the 

whole design 


process. Flow chart is the representation of transfer of control 
whereas the tree structure is the structure of subproblems 
generated in trying to achieve the goal. 

5,2 Tree Structure for the Design : 

The full tree structure for the present expert system is 
shown in Fig. 5.2 . Each goal (query) may give rise to several 
subgoals (subqueries). These subqueries are to be satisfied from 
left to right in the tree-structure. Each query or subquery 
corresponds to a predicate name in the program. For each predicate 
the rules present in the expert system have been indicated through 
their names in the elliptical boxes below each subquery or query . 
In tree structure 1 below each predicate three symbols appear t 
letter Q, letter F and A . The Q signifies a question-rule which 
causes the system to issue a question for user response. The F 
indicates that the fact exists for the predicate in 

. 7s " ^ - 



consideration. The presence of A signifies that the value of this 
predicate is already known and stored or the tree structure for 
this predicate is already described somewhere else in Fig. 5.2. To 
save the space, many times two or more rules for a predicate have 
been given in one tree diagram if the logic for them is the same. 

The tree structure shown in Fig. 5. 2 is a static 
template in which each rule and fact is shown only once. The 
actual tree of goals that will be generated during a computation 
will be dynamic; several instances of the same rule may occur in 
this tree, and the tree itself will vary from session to session 
depending on the user’s responses. However, the sequence of steps 
that it will perform will be the same as that shown in the 
flow chart of Fig. 5.1. 

Let us examine the initial part of the tree structure in 

detail : 

The top-level goal is ’complete’, for which two rules 
cl59 and cl60 exist. The selection of one of these rules depend 
upon the value of predicate ’ select-condenser ’ (whether 
vertical-upf low or vertical-downf low) . From top-level query in 
cl59, following sub-queries are generated: 

-select-condenser 

-cond-type 

-tube-diameter 

-flooding-mass-velocity 

-mass-velocity-tubes 

-ask-user 

-area-req. 

However, if the rule cl59 fails at the select condense^^ 
the rule cl60 is tried and the new goal ’area-req’ is generated. 
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Consider the rule cl59, here the current goal is 
select-condenser’ which has rules c8 to cll. It generates three 
further subqueries : 

-cond-type 

-t-or-p-condenser 

-ask-user 

For cond-type 8 rules {r9 to rl7) exist and they are 
tried in the order in which they appear in data base. Rule r9-rl0 
generates two subqueries 'pressure’ (rl) and 'retract-all’. rl 
being a question-rule, it issues a question to user and depending 
on user’s response it either succeeds or fails if ‘pressure’ 
succeeds, ‘retract-all’ predicate causes all other rules for 
‘cond-type’ to be unusuable in any further reference and the fact 
for ‘cond-type’ is the only thing which is accessible for later 
use. Rules rll-rl7 are self-explanatory on similar lines. 

‘t-or-p-condenser’ has again got a question rule and its 
third predicate is ask-user which causes issuing of question 
regarding type of condenser. Thus ‘select-condenser’ goal is made 

IT 

true if every thing goes in a proper way. 

The second subquery in rule cl59 is ‘cond-type’ which 

has already been made true while satisfying the goal 

‘select-condenser’. The next generated subgoal ‘tube diameter’ has 
got one question rule il and one fact i2. If the rule il fails 

then the ‘tube-diameter’ is satisfied from the fact i2. The other 

remaining subgoals are self-explanatory. The last subquery in rule 
(cl59) is ‘area-req’ which has three subgoals: 

‘area-req- condensing’ , ‘area-desup’ and ‘ area- subcool’ . 

The ' area-req-condensing’ generates a new goal 

* step-inlet -temp-area-upto-step’ , which makes use of recursion in 
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getting a true value. There are three rules c75, c76 and c77 for 
this predicate. The rule c75 is the boundary condition which 
causes recursion to end. The rule c75 generates a goal 
’coolant-temperature-at-the-start-of -condensation’ for which two 
rules c79 and c80 are present. 

The rules c76 and oil follow the same logic and the 
choice in them is made based upon the new subquery ’flow-type. 
The other subqueries are : 

-dely 

-step-inlet-temp-area-upto-st,ep (causes recursion) 

-heat - load-per- step 

-fluid-name 

-mass-f low-rate-cold 

-specific-heat 

- over al 1 -htc -condens ing 

-condensing-temp 

-log-mtd 

And thus the tree for the top goal ’complete’ keeps on 
growing and it can be understood based on the features explained 
till now. 

5.3 Structuring of the Knowledge Base : 

5.3.1 Type of Data Structure : 

The structure of the rules to a great extent depends 
upon the type of data structure used. For example, consider the 
fluid properties ( property- 1, property- 2. ... etc) for the fluids 
fluid-l and fluid-2. Now there can be various types of data 
structure : 
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(I) 

(f luid-1 

property -1 

<vai> ) 


(fluid-2 

property-2 

<val> ) 


( property -1 

fluid-1 

<val> ) 

(II) 

( property- 1 

fluid-2 

<val> ) 


(property- 2 

fluid-1 

<val> ) 


(fluid-1 property-1 .. .property-n <val l>...<val n>) 
(fluid-2 property-1 .. .property-n <val l>...<valn>) 


(property-l-and-2 fluid-1 

<val 

1> 

<val 

2>) 

(property-l-and-2 fluid-2 

<val 

1> 

<val 

2>) 

(property-l-of -fluid-1 

<val 

11>) 



(property-l-of -fluid-2 

<val 

12>) 



(property-2-of -fluid-1 

<val 

21 >) 




and so on. 

Let us consider the types (I) and (II). Type I is 
centered around fluids (fluid-1, fluid-2 etc) while type-II is 
centered around the properties (density, temperature etc here 
written as property-1, property -2 etc). Since the number of 
fluids is small in our case, viz., 2, we get just two predicates 
and so the number of explanations, questions to be provided etc 
gets reduced in number where as type II providing a large number 
of predicates gives flexibility of association of different 
questions, explanations etc with each one of them. The problem of 
type I exists also in type III though it saves space. 

It is important to understand that there has to be a 
practical limit beyond which additional flexibility becomes a 
burden rather than a help. For example, in type V the 
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expianation for (say) density would have to be repeated over and 
over again. Wherever the same explanation is sufficient for 
several things, they can be combined in one predicate as for 
example in type IV. 

The type IV is useful when the fluid is undergoing 
phase change and we need the properties of both the phases . Then 
the <val 1> can be property of liquid phase and <val 2> of vapour 
phase. In the present expert system both type II and type IV have 
been used. 

5.3.2 Different Predicates Used : 

The computational predicates take their arguments and 
return true (T) or false (Nil) according to the condition e.g.. 
<, >,le, ge etc. (le = less than or equal to and ge = greater than 
or equal to). Most of the predicates defined for this system take 

one argument which is a numerical value e.g. heat-rate, 

■ % 

no-of -tubes, shell-diameter. Such predicates essentially have a 
rule, and most of the time they are inferred from a rule. 

Other types of predicates are those which have more than 
one argument. The need for defining such predicates arises because 
of change in references and requirement of storing many 

properties, e.g. : 

♦ 

(specific-heat fluid-1 <temp-val> <val 1>) 
(specific-heat fluid-2 <temp-val> <val 2>) 

The built-in predicates in VIDHI like ’ask-user’ and 
’ retract-all ’ have their own significance and have been explained 
■ earlier. ■ 
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5.3.3 Types of Rules : 


The rules used here are mos'tly of "the type which makes 
use of an expression for assigning value to a variable argument of 
the predicate : 

A rule of the form 

(defasrt sll9 ( subcool-ht-load ?qs) <- 
(mass-flow-rate condensation-side ?mh) 

(fluid-name hot-fluid ?name-h) 

(outlet-temp ?name-h ?tho) 

(condensing-temp ?name-h ?tsat) 

(specific-heat ?name-h ?cpl ?cpg) 

(= ?<js (product ?mh ?cpl (diff ?tsat ?tho)))) 

is equivalent to the expression 

9 1 . =m, X c , X (T - T ) 

sub >1 pi 'ho sol ' 

This can be asserted as a logic expression in other 
rules as follows; 

(defasrt S63 (area-subcool ?acool) <- 
(subcool-ht-load ?qs) 

(overall-htc-subcool ?uo) 

(= ?acool (quotient ?qs (product ?uo...)))) 

The second type of rule is one which takes purely 
logical decisions, consider the rule; 

(defasrt c3 (shell-side-fluid ?name-s) <- 
(cond type shell-side-condensation) 

(condensing-fluid ?name-s) 

(retract-all shell-side-fluid)) 

(defasrt c4 (shell-side-f luid ?name-s)<- 
(cond-type tube-side-condensation) 
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(fluid-name cold-fluid ?name-s) 

(retract-all shell-side-fluid)) 

The first rule says that if shell-side-condensation is 
going on then the condensing-f luid is the shell-side-fluid and the 
second rule says that if it is case of tube-side-condensation then 
the .cold-fluid (i.e. non-condensing fluid) is the 
shell-side-fluid. 

It is to be noted that in the first type of rules the 
ordering of subqueries in antecedent is important , because all the 
variables used in a computational functions or predicates should 
have bindings by the previous subqueries, 

5 . 4 Encoding ^ New Rule ; 

Let us consider a brief example illustrating how the 
available knowledge is converted into a rule. Consider the 
following knowledge available about pitch-layout : 

1 . The most commonly used pitch-layout is the 
staggered-30. 

2. If the restriction on shell-side pressure-drop is 
strict, then inline-90 is the suitable type. 

3. If the requirement of mechanical cleaning is 
frequent then staggered-45 is the most suitable 
type. 

The first step is to arrange these statements in the 
order of importance. 

Along with this, the important parameters are selected 
as predicates. So we have : 
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Order of Importance Serial -number 

1 2 

2 3 


3 


1 


Predicates 

( pres sure- drop- strict ) 

(mechanical-cleaning- 
f requent ) 


For the type of pitch-layout we can have predicate 
’ pitch- layout ’ . 

Now using the VIDHI syntax, rules for ’ pitch- layout ’ can 
be written as follows : 

(defasrt f44 (mechanical-cleaning-frequent) <- 
(ask-user (source) (target) 

(question Is the requirement of shell side 
mechanical -cleaning-f requent ) ) ) 

(defasrt f 45 (pressure-drop-strict) <- 
(ask-user (source) (target) 

(question Is the limitation on shell side pressure-drop 
strict))) 

(defasrt f46 (pitch-layout ?x) <- 
(pressure-drop-strict ) 

(= ?x inline-90) 

(rertract-all pitch- layout ) ) 

(defasrt f47 (pitch-layout ?x) <- 
( mechanical -cleaning-f requent ) 

(= ?x staggered-45) 

(retract-all pitch-layout)) 

(defasrt f 48 (pitch-layout ?x) <- 
{= ?x staggered -30) ) 
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It can be noted that the rules f46 and f47 cannot be 
interchanged since if the user wants to satisfy the requirements 
of both the Predicates (pressure-drop-strict) and 
(mechanical-cleaning-frequent), then we want f46 to get control 


rather than f47. 
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6- RESULTS AND DISCUSSION 


6 . 1 Hqw tQ Use the System : 

If the decision regarding the use of shell-and-tube 
condenser having an E-type shell has been taken, the present 
system can be used to design it. The system is very interactive 
and user friendly. It issues the questions to the user at every 
step giving sufficient help wherever possible. 

The questions asked by the system are of two types : 
a. Mandatory Questions and b. Optional Questions 

Mandatory questions collect that vital information from 
the user without which nobody can even think of designing a 
condenser. If the user responds incorrectly to this information 
or, is not able to provide it, the system can’t proceed to design 
the condenser. The mandatory information needed has been listed in 
Table 6.1 and must be collected before proceeding to execution. 

Optional questions are not so critical and even if user 
responds to these questions saying "dontknow”, still the design 
can proceed. These optional questions are issued to the user to 
enable him getting a tailor-made design according to his need. For 
example, user may prefer only one particular shell-diameter or may 
have an upper limit on its value; in that case he should specify 
it. Also, if the user before hand .has some inf onnation regarding 
design, (in the second iteration he definitely has some) e.g., 
number of tubes, number of baffles, etc. , he should specify it. 
This will save unnecessary time used in computing these 
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TABLE 6 . 1 

Questions 


Question 


Condensing fluid Coolant 


Name of the fluid 


Specific heat 

J/kg-K 

<value-liquid><value-gas> 

<value> 

Thermal conductivity 

W/m-K 

<value-liquid><value-gas> 

<value> 

Thermal conductivity of 

tube material 

Inlet temperature 

W/m-K 

°C 

<value> 

<value> 

Outlet temperature 

"c 

<value> 

<value> 

Condensing temperature 

°C 

<value> 

<value> 

Mass flow rate 

kg/s 

<value> 

<value> 

Absolute viscosity at 

average temperature 

cP 


<value> 

(Liquids only) 

Absolute viscosity at 

some other temperature 

cP 


<value> 

Absolute viscosity at 

condensing temperature 

cP 

<value-liquid> <value-gas> 


Density 

kg/m^ 

<value-liquid><value-gas> 


Condensing pressure 

kPa 

< value > 


Enthalpy of vaporization J/kg 

<value> 



Approximate overall heat- 

transfer coefficient W/m*-'*C 
Viscosity of condensate cP 
Aspect ratio 

Flow-direction 
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TABLE 6 . 2 


Optional Quest, ions 


Question 

Shell diameter 

mm 

<value>/dontknow 

Tube diameter 

mm 

<value>/dontknow 

Area of heat transfer 

2 

m 

<value>/dontknow 

Number of baffles 


<value> / dontknow 
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6 . 2 §.asml e 


Two sample sessions are presented in the Appendij^o 
illustrate the behaviour of expert system under different type of 
responses provided by the user. 

Session 1 illustrates the case where a mandatory 
question is not answered. The system is not able to proceed beyond 
that point. 

Session 2 illustrates the normal behaviour of the system 
as expected. The user answers the questions in a proper way after 
seeking the explanations whenever he is in doubt. 

6 • 3 Limitations of the Present System and Scone for Future Work 

The limitations of the present system are because of two 
reasons. One set of limitations exist because the VIDHI itself is 
undergoing development, which of course requires further work on 
the part of a Computer Scientist. The other set of limitations is 
due to insufficient knowledge base which still requires additional 
data and literature survey. These can be summarised as follows ; 

a. The rules for the calculation of pressure drop in two-phase 
zone are not included. So the pressure rating of the design is not 
a part of this thesis. 

b. Present work considers only the case of single shell-pass and 
single tube-pass condensers. For the case of multi tube -pass, the 
step-wise method in two-phase zone will have to be modified 
considerably and becomes very complicated. However , using an 
average condensation coefficient, this extension can be easily 
incorporated. 
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c. The present system is only for E-type shell. 

d. The present system lacks the provision of changing an already 
given value of a particular design parameter. This arises because 
of inability of declaring data dependencies in the ’Shell’ and 
absence- of such provision in the tree structure. 

All the above limitations provide the scope for future 

work. 

6.4 Conclusions : 

The expert system developed here makes use of the ’Shell 
- VIDHI’ (Sangal, 1988) . It is based on logic programming and is 
highly interactive. It is capable of designing a vertical 
condenser from heat transfer point of view. The rules involving 
the shell side computations make use of ’ Bell -Del aware Method’ in 
the rating program. In the condensation zone, the stepwise or 
successive summation method is used to take into account the 
variation in condensation coefficient due to the change in vapour 
fraction of the mixture. After the rules for the two-phase 
pressure drop calculations are added, the system will be able to 
completely design a single tube, single shell pass vertical 
condenser using E-type of shell from the thermal performance point 
of view. 

Additional rules can be easily added to include other 
types of multipass shell -and-tube condensers. After some field 
testing, the system can be put into commercial use. Thus, the 
present work marks only a beginning towards this ultimate goal. 
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APPEKDIX 


SAMPLE SESSION 1 


% ~sangal/vidhi .exe 

14. load cini 
Cload cinil 
Cload compD 
Cload listU 
Cload list13 
Cload cond3 
Cload alv13 
Cload cnds13 
Cload cnds23 
Cload cnds33 
Cload cnds43 , 

Cload- cnds53 
Cload cnds63 
Cload itercond3 
Cload fixedcond3 
Cload vidhiupdate . 13 
t 

1S. design 

("what is the pressure of condensing vapour * 
high , very-low, or intermediate" > intermediate 
("Is the vapour very-corrosive , corrosive , or non-corrosive" ) corrosive 
(is the condensation total or part ial ) total 
("We have a case of" tube-side-condensation " 

*Now please select one of the types after 
studying the following tabled 



horizontal 

vertical-downflow 

^ 

vertical-upf low 

single-component 

g 

av 

9 

av 

f 

b 

vapour 







multi-component 

f 

b 

g 

b 

f 

C 

vapour 







subcooled- 

P 

K 

9 

av 

K 

K 

condensate 







pressure-drop 







high 

g 

b 

9 

c 

K 

C 

low 

P 

b 

f 

'C'. ■ 

g 

C 

coolant 







liquid 

g 


g 


9 


gas 

9 


9 


9 


boiling 

9 


9 


9 


acceptabi lityi g® 

good , 

pspoor r 

X® not 

acceptable or 

not 

recommended 


f»fair 
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predictability* av=average 257. 

b =fair *^<50% 
c =poor *^>50% 

X =no method or not recommended 

horizontal condensers have possibility of slugging, 
vertical downflow handles dirty or polymerizing vapours, 
please answer on right hand side") vertical— upf low 
(Give the value of the tube outside diameter in mm)what 


( the recomended tube dimensions are *- 


tube o.d 

mm 

wall thickness 
mm 

tube i.d 

mm 

outside surface 
sq m / m 

6.0 

0.5 

5.0 

0.019 

8.0 

1.5 

5.0 

0.025 

10.0 

1.5 

7.0 

0.031 

14.0 

2.0 

10.0 

0.044 

18.0 

2.0 

14.0 

0.057 

20.0 

2.0 

16.0 

0.063 

25.0 

2.5 

20.0 

0.079 

30.0 

2.5 

25.0 

0.094 

38.0 

2.5 

33.0 

0.119 

44.5 

2.5 

39.5 

0.139 

51.0 

2.5 

46.0 

0. 160 


select a suitable tube diameter or tube wall thickness from the above 
table . generally, tubes below 20.0 mm o.d are not recommended since 
they have a large pressure drop, smaller diameter tubes give larger 
heat transfer area for the same shell diameter, 
you may say dontknow also. )>38.0 
(what is the name of condensing fluidlsteam 
("what is the condensing temperature of" steam " 
in degree Celsius at condensing pressure")what 

(This is a mandatory question) >100.0 
("Give the absolute viscosity of liquid and vapour 
phase respectively of" steam "in centipoise at 
temperature" 100.0 deg celsius)0. 2^8 0.0126 
("please specify the liquid and vapour phase density 
of" steam at temperature 100.0 in kg/m3)958.3 0.597 
("what is the value of surface-tension coefficient 
for liquid-condensate in mN/m, or , dyne/cm")58.0 
(what is the taper of tube-end in degrees from horizontal)60.0^^^^^^^^^ ^^^^^^^^^ 
(Give the value of the tube-wall-thickness in mmlwhat 

( the recomended tube dimensions are *- 
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tube o.d 

mm 


wall thickness tube i.d 
mm mm 


outside surface 
m / m 


0.5 5.0 

1.5 5.0 

10.0 1.5 7,0 

1^*0 2.0 10.0 

18.0 2.0 14.0 

20.0 2.0 16.0 

25.0 2.5 20.0 

20.0 2.5 25.0 

28«0 2.5 33.0 

2.5 39.5 

51.0 2.5 46.0 


0.019 
0.025 
0.031 
0.044 
0.057 
0.063 
0.079 
0.094 
0. 119 
0. 139 
0 • 1 60 


select a suitable tube diameter or tube wall thiru 

table . generally, tubes below 20.0 mm o.d are not above 

they brave a lar-ge- -prese-ure- -sma-Mer-djetmat-er ^ '"®^°"''’>ended since 

heat transfer area for the same shell diameter. — 

you may say dontknow also. )>2.5 
(What is the name of the cold—f luidlwater 
(what is the mass-flow— rate of the condensing f lui ^ • 

(Give the value of the shell diameter in mmldontknr, *^’9''s>3.5 
(Give the area of heat transfer in sq m)what 


(if this the first round of calculations or 
if you know the value of overall heat transfer coi»<‘x‘- • 

W/ sq m K, say *dontknow* else give the value from 
initial round of calculation )>ontknow 
("Give the value of overall heat transfer coefficio ♦ • 
W/sq-m-K from previous experience or from 
explanation available"}dontknow 
(what is the taper of tube— end in degrees from hor-i 
("what is the flow-directions of the streams: 2ontal)60.0 

CO— current, counter") counter 


("How many number of steps should be choosen to 
determine required condensing area by step-wise 
(" what is the value of enthalpy of vaporization "'®‘'bod")4 
of condensing fluid in J /kg" >2256700.0 


('please specify the specific heats of liquid and va 

phase respectively of" steam in j/ka-u?^oo/^ 
(Please specify the inlet-temp of steam in deg-celsi^ 

(Please specify the outlet-temperature of steam In 110.0 
(Please specify the inlet-temp of water in deg-celsi^^jop 
(Please specify the outlet-temperature of water in . 

(please specify the specific-heat of water in j/kg k>4i7A 
(please specify the conductivity of tube material 
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this is a table for comparision of the results 


the heat transfer area required sq m = (♦var areql 

the actual heat transfer area available sq m = ($var a-o> 

the desuperheating area required sq n» = ($var a-desup> 

the condensation area required sq m = ($var aeon) 

the subcooling area required sq m = ($var acool) 

shell side fluid « ($var name-s) 

tube side fluid *= steam 

no of baffles = (*var n-b) 

no of tubes »= (♦var n-t) 

tube side passes *= 1 

shell side passes = ($var ns) 


the different correction-factors for heat transfer are 


window correction factor = ($var j-c) 

baffle leakage correction factor = (♦var j— 1) 
bundle bypass correction factor = <$-var j-b) 

adverse temperature gradient correction factor = <$var J-r) 

total correction factor = Error! Non-number to add (*var j-c) 
<1>! (reset) 

16. exit 

■/. 


NOTE! Here the User did not specify either Approximate Area 
or approximate Overall Heat Transfer Coefficient, so 
the Design could not be completed. 
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SAHPLE SESSION 2 


7 . '^sangal/vidhi .exe 

14. load cini 
Cload cini3 
Cload comp3 
Cload li'st3 
Cload listl3 
Cload cond3 
Cload alv13 
Cload cnds13 
Cload cr»ds23 
Cload cnd.s33 
Cload cnds4Il 
Cload cndsI53 
Cload cnds63 
Cload itercondH 
Cload fixedcondH 
Cload vidhiupdate.il] 
t 

15. design 

("what is the pressure of condensing vapour s 
high, very-low, or intermediate" ) intermediate 
("Is the vapour very-cor ros i ve , cor rosi ve , or non-corrosive") corrosive 
(is the condensation total or partial ) total 
("We have a case of" tube-side-condensation " 

*Now please select one of the types aftpr 
studying the following table* 



horizontal 

vertical-downflow 

■* — -- - — * 

vertical-upf low 

single-component 

vapour 

g 

av 

9 

av 

f 

b 

multi— component 

f 

b 

9 

b 

f 

, c 

vapour 







subcooled- 

P 

X 

9 

av 

K 

K 

condensate 







pressure-drop 







high 

9 

b 

9 

■ c 

X 


low 

P 

b 

f 

, c 

9 

c 

coolant 







liquid 

9 


9 


9 


gas 

9 


9 


9 


boiling 

9 


9 


9 
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acceptability! g=good , p=poor , k'= not acceptable or not recommended 

f !=fai r 

predictability* av’==average 25% 

b “fair '^<50% 
c “poor *^>50% 

X =no method or not recommended 

horizontal condensers have possibility of slugging, 
vertical downflow handles dirty or polymerizing vapours, 
please, answer on right hand side‘'>vertical-downf low 
("what is the flow-directions of the streams* 
co-current, counter''ywhat 

(Of course the system is developed for cross 

-flow E-type of shells, so if the both fluids enter on the same-side 
it is a case of co-current, else counter. 

This is a mandatory question) >counter 

("How many number of steps should be choosen to 

determine required condensing area by step-wise method")4 
(what is the name of condensing fluid>what 

(this is a mandatory question) >steam 

(what is the mass-flow-rate of the condensing fluid in kg/s)3.5 
(" what is the value of enthalpy of vaporization 
of condensing fluid in J/kg")2256700.0 
(What is the name of the cold-f luid)water 
("please specify the specific heats of liquid and vapour 

phase respectively of" steam in j /kg-K )4E20.0 2030.0 
(Please specify the inlet-temp of steam in deg-celsius)110.0 
(Please specify the outlet-temperature of steam in deg celsius)90.0 
("what is the condensing temperature of" steam " 

in degree Celsius at condensing pressure" ) 100.0 
(Please specify the inlet-temp of water in deg-celsius )32.0 
(Please specify the outlet-temperature of water in deg celsius)40.0 
(please specify the specific-heat of water in j/kg k)4178.5 
(please specify the conductivity of tube material in W/m- K)38.0 
(Give the value of the shell diameter in mm)what 

( if there is a restriction on 
the shell-diameter give the value else say 
dontknow) >dontknow 

(Give the area of heat transfer in sq m)what 

(if this the first round of calculations or 
if you know the value of overall heat transfer coefficient in 
W/ sq m K, say ’dontknow* else give the value from the 
initial round of calculation )>dontknow 
("Give the value of overall heat transfer coefficient in 
W/sq-m-K from previous experience or from 
explanation available")what 

(some typical ranges are given in 
the following table — 
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¥apo'iir 

coolant 

U,W/sq-m-K 

alcohol 

water 

50-1100 

dowtherm 

tall oil 

40-450 

dowtherm 

dowtherm 

450-680 

high-biol ing-hydro carbons 

under vacuum 

water 

100-280 

low-biol ing-hydrocarbons 

water 

450-1140 

organic solvents 

water 

550-1140 

hydrocarbons 

oil 

140-230 

kerosene 

water 

170-370 

kerosene 

oil 

110-170 

naphtha 

water 

280-430 

naphtha 

oil 

110-170 

steam 

feed water 

1800-5700 

vegetable-oils 

water 

110-280 

organic steam, azeotrope 

water 

220-450 

air-coolers 

steam 

ai r 

730-800 

— - - - * 

ai-r- - - — ‘ 

_550-680- — 

light-hydrocarbon 

air 

450-540 

light-naphtha 

air 

400-450 

f reons 

air 

340-450 

heavy-naphtha 

ai r 

340-400 


please answer on right hand side 

This is a mandatory quest ion ) >1800.0 

<Is the limitation on shell side pressure-drop-strict)n 

<Is the requirement of shell side mechanical-cleaning-f requent)n 

(Give the value of the tube outside diameter in romlwhat 

< the recomended tube dimensions are s- 

tube o.d 

wall thickness tube i.d 

outside surface 

IHflU 

mm 

mm 

sq m / m 

. 

6.0 

0.5 

5.0 

0.019 

8.0 

1.5 

5.0 

0.025 

10.0 

1.5 

7.0 

0.031 

14.0 

2.0 

10.0 

0.044 

18.0 

2.0 

14.0 

0.057 

20.0 

2.0 

16.0 

0.063 

25.0 

2.5 

20.0 

0.079 

30.0 

2.5 

25.0 

0.094 

38.0 

2.5 

33.0 

0.119 

44.5 

2.5 

39.5 

0.139 

51.0 

2.5 

46.0 

0.160 
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select a suitable tube diameter or tube wall thickness from the above 
table . generally, tubes below 20.0 mm o.d are not recommended since 
they have a large pressure drop, smaller diameter tubes give larger 
heat transfer area for the same shell diameter, 
you may say dontknow also. >>3S.O 
(Give the value of the aspect ratioM2.0 
(Is the steam foulingln 

(Give the inlet-pressure of the steam in kpa)101.3 
(Is a si ight ly-greater-diameter of the shell tolerablely 
(Give the inlet-pressure of the water in kpa)101.3 
(Give the number of baffles )what 


( the no of baffles decide the shell side 
reynolds number, lower value gives a lower Reynolds Number 
and lower pressure drop. 

if the value of baffle leakage correction factor is less than 
0.7 reduce the number of baffles. This increases the cross flow 
area at shell centre and hence the correction factor. )>dontknow 
(“Give the absolute viscosity of" water " 
at" 36.0 deg celsius in cP)0.710 
(Give the thermal conductivity of the 


(please classify the water "in any of 
following categories s- 
(Dlight-liquid 

(2) medium-l iquid 

(3) heavy-hot-1 i quid 

(4) heavy-cold— 1 iquid 

(5) very-heavy-hot- 1 iquid 

(6) very-heavy— cold-1 iquid 
(7>gas-at-a-pressure < 200 kpa 

(8) gas-at-a-pressure < 1000 kpa 

(9) gas-at-a-pressure < 10000 kpa 


water 

the 


in W / m K)0.626 


light 

medium 

heavy-hot 

heavy-cold 

very-heavy- 

very-heavy- 

gas-200 

gas- 103 

gas- 104 


hot 

cold 


the gas-pressure is the absolute pressure 

the liquid has to be classified on the basis of viscosity 
answer by typing the word on right hand side? this is a 
mandatory question")light 

("Give the absolute viscosity of the" water "at some other 
temperature. Give the temperature in deg celsius followed 
by the viscosity in cP.")95.0 0.299 
(Give the value of the tube-wall-thickness in mm)2.5 . 

("Give the absolute viscosity of liquid and vapour 
phase respectively of" steam "in centipoise at 
temperature" 100.0 deg celsius)0.278 0.0126 
("please specify the conductivity of liquid and 
vapour phase respectively 
of" steam in W/m-K 10.679 0.025 
("please specify the liquid and vapour phase density 
of" steam at temperature 100.0 in kg/m3)958.3 0.597 
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(please classify the steairi "in any of 
following categories s- 

(1)1 ight-1 iquid 
(Elmedium-liquid 

(3) heavy-hot-1 iquid 

( 4 ) heavy-cold-1 iquid 

(5) very-heavy-hot-l iquid 

(6) very-hea vy-co Id-liquid 

(7) gas-at-a-pressure < 200 kpa 
<8)gas-at-a-pressure < 1000 kpa 
(9) gas-at-a-pressure < 10000 kpa 


the 

light 

medium 

heavy-hot 

heavy-cold 

very-heavy-hot 

very-heavy-cold 

gas -200 

gas-103 

gas-104 


the gas-pressure is the absolute pressure 

the liquid has to be classified on the basis of viscosity 
answer by typing the word on right hand side; this is a 
mandatory quest ion" ) gas-200 
("please classify the condensed fluid 
of the following categories*-* 

(1) light-liquid 

(2) medium-liquid 

(3) heavy-hot-l iquid 

(4) heavy— cold-1 iquid 

(5) very-heavy— hot-1 iquid 

(6) very-heavy-cold-liquid 


in any one 

light 

medium 

heavy— hot 

heavy-cold 

ve ry-heavy— hot 

very-heavy-cold 


condensed liquid has to be classified on the basis of 
viscosity , answer on right hand side")light 
("what is the value of viscosity of condensate ( 1 iquid) 
in centipoise at temperature" 95.0 deg celsius)0.299 
("please classify the vapour in any one of the following 
categories!-* 

C11 vapour at a pressure < 200 kpa gas-200 

C2I1 vapour at a pressure < 1000 kpa gas-103 

C33 vapour at a pressure < 10000 kpa gas-104 

the gas pressure is absolute pressure 
answer on right hand side")gas-200 
("what is the value of viscosity of condensing vapour 
in centipoise at temperature" 105.0)0.0125 


yes 

Should I try for another answer(y/n) !n 


this is a table for comparision of the results 


the heat transfer area required sq m = 68.74666349468933 
the actual heat transfer area available sq m = 70.6693321090584 
the desuperheating area required sq m = 5.345516950502597 
the condensation area required sq m = 61.86690341596033 
the subcooling area required sq m = 1.534243128226397 
shell side fluid = water 
^ - t4ibe 5-i^ 44uid 
no of baffles =10 
no of tubes =63 
tube side passes =1 
shell side passes =1 
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window correction factor =. n coDc:-rr=n«^^Iir“ “ 


i>aff le leakage cor ection fact ^'■S^SSysSOSSOSSasi 

bundle bypael correcli^rflrt^r : 
adverse temperature gradient correction ff r^f 
ai correction factor ^ O.4415S328520976O3 


the shell side parameters 


shell diameter mm = 495.0370S394 oap7i 

length of the shell mm = 5940 . 445007 SSa 3 AA 
no of baffles = 10 - ^4:>oo728S366 

central baffle spacing mm = 53 1 . 03978 i 854 Aooe; 
percentage of baffle cut = 36 . 363 A 3 A 3977 P 7 pa ^ 
centriangle of baffle cut deg ^ 
upper centriangle of baffle cut 

bundle to shell inside diametej clelrencrm^- 

aspect ratio = is.o learence mm - 96.08870522807266 

cross flow area at the she! 

Tr^cirion no of tubes in pure crnc:c - a ^ 

bundle to shell byoass area so Ja™ 0- 3870497349629656 

shell to baffle leajlje irJa L " ' '>0937.16920777183 
tube to baffle holftelkS^ arf ' 2322.487039472815 

nS "V 5.080148235762788 

no ot sealing strip pairs =0 

no of shell side passes =1 


the tube side parameters 


name of the fluid = steam 

tube diameter mm = 38.0 

tube wall thickness mm == 2.5 

tube length mm = 5841.437600400968 

no of tubes = 63 / 

no of tube side passes = 1 

^tube bundle type = pull-through-f loating-type 

16. exit 
% 
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